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The interlayer exchange coupling mediated by valence band electrons in all-semiconductor IV- VI 
magnetic/nonmagnetic superlattices is studied theoretically. A 3D tight-binding model, accounting 
for the band and magnetic structure of the constituent superlattice components is used to calculate 
the spin-dependent part of the total electronic energy. The antiferromagnetic coupling between 
ferromagnetic layers in EuS/PbS superlattices is obtained, in agreement with the experimental 
evidences. The results obtained for the coupling between antiferromagnetic layers in EuTe/PbTe 
superlattices are also presented. 



I. INTRODUCTION 

Interlayer exchange coupling (lElC) was discovered in 
late 80-ties in Fe/Cr/Fe trilayersEI. Since then it has 
been observed in a variety of multilayer structures com- 
posed of alternating magnetic and nonmagnetic lay- 
ers. These studies concentrated on the coupling be- 
tween ferromagnetic, metallic layers, although separated 
by both metallic and insulating spacers. Thus, the fer- 
romagnetic character of the magnetic layers and the 
fact that in these structures the Fermi level was situ- 
ated in the region of high density of electronic states 
inhered to the theoretical models, which were designed 
to explain the origins of the TEC phenomena (seeEl and 
the references therein). Surprisingly enough, the lEC 
was also discovered in all-semiconductor superlattices 
(SLs). Moreover, the semiconductor SLs in which iJ. 
was first observed were the MnTe/CdTcfl, MnTe/ZnTefil 
and EuTe/PbTeO, all with antiferromagnetic layers. Re- 
cently, such coupling was also identified in semicon-, 
ductor multilayer strujctures with ferromagnetic EuSJj'lI 
and Ga(Mn)As layersB While the lEC in the trilayer 
Ga(Mn)As/Ga(Al)As/Ga(Mn)As with the high concen- 
tration of free carriersjj, can be, at least qualitatively, 
explained in terms of the models tailored for metallic 
systems, the other semiconductor structures exhibiting 
lEC call for a different approach. 

Several attempts to explain lEC in all-semiconductor 
structures have already been reported in the literature. 
Two models, in which the interlayer coupling is medi- 
ated by carriers localized on shallow impunties in the 
spacer region, were proposed for II- VI SLsp't2l. These 
models do not apply to IV-VI structures, with the 
PbTe and PbS spacers, since in lead chalcogenides la- 
calized shallow impurity states were never detectedtil. 
For these SLs, mechanisms of interlayer spin-spin in- 
teractions mediated by valence-band electrons were sug- 
gested. The calculations of the difference between to- 
tal electronic energies obtained for two different (Bl 
and B2 in Fig. ^ spin configurations of the SL, per- 
formed within a frame of a very simple ID tight-binding 



model, put first in evidence the significant role of the va- 
lence band electrons in lEC in all-aepiiconductor mag- 
netic/nonmagnetic layer structuresta. This role was 
further demonstrated for EuTe/PbTe/EuTc trilayers by 
Wilczyhski and Swirkowicz inO, where a 3D tight bind- 
ing model, still oversimplifying the band structure, was 
used. A different approach to the magnetic interlayer 
interactions mediated by jzalence electrons have been 
chosen by Dugaev et al,EJ. These authors studied 
the Blombergen-Rowland mechanism within the effec- 
tive mass approximation and obtained a ferromagnetic 
coupling between two magnetic impurities situated at 
the opposite interfaces of a narrow-gap IV-VI semicon- 
ductor spacer. As the experimentalW observed lEC in 
EuS/PbS SLs is antiferromagneticpQ, this means that 
the Blombergen-Rowland interactions are not dominat- 
ing lEC in these SLs. Still, because of the low concentra- 
tion of free carriers and the absence of shallow impurity 
states in PbS, the lEC in these SLs is most likely me- 
diated by valence electrons. In this situation, the total 
energy calculations, which do not focus on a particular 
interaction mechanism but account globally for the spin- 
dependent structure of valence bands, seem to consti- 
tute the most approjjtiate approach. The calculations of 
this type, reported inllj andE3, were clearly oversimplified 
and performed for a different spin structure than that of 
EuS/PbS (001) SL. 



In this paper we present the results of refined 3D total 
energy calculations, which take into account the crystal 
and the band structures of the SLs' component materials. 
The tight-binding model with its assumptions and the 
results for three different spin structures, corresponding 
to all experimentally studied IV-VI semiconductor mag- 
netic/nonmagnetic SLs, are presented in Section III. In 
Section II the magnetic and band structures of the con- 
stituent materials are described and the SL geometry is 
specified. The comparison with experimental data and 
the conclusions are presented in the last Section. 
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II. CONSTITUENT MATERIALS AND THE 
SUPERLATTICES GEOMETRY 



All the components of the EuX/PbX SLs, where X=Te 
or S, crystallize in the rock-salt structure. Bulk PbS and 
PbTe are narrow gap nonmagnetic semiconductors with 
very similar band structures. They both have a direct 
energy gap between the p-anion valence band and the 
p-cation conduction band at the point L of the BZEj. 
It is well known that in both these lead chalcogenides 
the spin-orbit terms are important for the detailed de- 
scriptioji of the energy bands in the vicinity of the L 
pointjlj. One notices, however, that the spin-orbit cor- 
rections affect predominantly the conduction bands orig- 
inating from the states of the heavy Pb atoms. 

Bulk EuS is a classical Heisenberg ferromagnet with 
the Curie temperature 16.6 KLD. Bulk EuTe exhibits 
type II AFM structure with the Neel temperature 9.6 
KEZI. In EuTe the spins of Eu ions are ferromagnetically 
ordered in (lll)-type planes, which in turn are coupled 
antiferromagnetically to one another. All the Eu chalco- 
genides are semi-insulating, large gap semiconductors. 
The results of the nonrelativistic APW., calculations of 
the EuS spin-polarized band structureta show the nar- 
row f(t) bands situated in the energy gap between the 
valence band, formed essentially of anion p states, and 
the conduction band, built mostly of cation d states. The 
valence band maximum is situated at the center of the 
Brillouin zone (BZ) and the conduction band minimum 
at the point X. The spin splitting of the valence band re- 
sults predominantly from the spin dependent mixing of 
p anion and f cation states, whereas that of the conduc- 
tion band is mostly due to f-d and s-d on-site direct ex- 
change. Much less is known about the EuTe band struc- 
ture. The optical experiments performed at T=300 K 
indicatdlj, that in the paramagnetic phase of EuTe the 
f-d gap is somewhat larger than in EuS, in agreement 
with the general trends visible in the experimental datacll 
and the results, of APW calculations for other europium 
chalcogenidesEa. These trends seem not to be followed in 
the recentEj calculations of EuTe band structure, focused 
predominantly on the conduction bands. 

Two types of EuS/PbS SLs were experimentally stud- 
ied, one grown on KCl substrate along [001] and the other 
on BaF2 along the [111] crystallographic axis. The mea- 
surements show that in both cases the ferromagnetically 
ordered Eu spins within each magneticp4ayer lie in the 
planes perpendicular to the growth axiaHJ. In the (001) 
structures each atomic monolayer consists of both anions 
and cations, with the monolayers a/2 apart, where a is 
the cubic lattice constant. The schematic view of two 
such monolayers is presented in Fig. |^. The distances 
between the cation and its four in-plane nearest neigh- 
bors (anions) and four in-plane next nearest neighbors 
(cations) are also shown in the figure. The spin struc- 
ture of the magnetic (001) SL, which corresponds to the 
observed antiferromagnetic interlayer coupling is shown 
schematically in Fig. |lj A2. 
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FIG. 1: The correlated co-linear spin structures for (A) fer- 
romagnetic and (B) antiferromagnetic layers. For the "in- 
phase" spin structures Al and Bl the magnetic period is equal 
to the chemical one, for the "out-of-phase" spin structures A2 
and B2 the magnetic period is two times larger. 



In the case of (111) EuS/PbS SLs lEC has not been 
yet observed. These SLs have the same crystallographic 
structure as the experimentally studied EuTe/PbTe SLs 
grown along the [111] axis. With this growth direc- 
tion, the subsequent (111) atomic monolayers are built 
either of anions or of cations, in alternation. The dis- 
tance between cation and anion monolayers is aV3/6. 
A schematic view of three successive cation layers is pre- 
sented in Fig. ||- the analogous anion sublattice is shifted 
by a^/3/2 along the [111] direction and is not shown in 
the figure for the clarity reasons. 

In EuTc/PbTe SLs the neutronographic measurements 
show that the AFM Type II structure is preserved in each 
EuTe layer, but the FM spin sheets form exclusively on 
the (111) planes parallel to the layers. This is in contrast 
to bulk EuTe, where there are four symmetry-equivalent 
Type II AFM arrangements in which the ferromagnetic 
spin sheets form on the {111}, {HI}, {HI}: or {111} 
plane families. Moreover, for the nonmagnetic spacers 
thin enough, the satellites observed in neutronographic 
spectra clearly indicate the existence of some long range 
order proving that the spins in consecutive magnetic lay- 
ers are not randomly oriented, but_tend to align along the 
same direction in a correlated wayd. Although for the an- 
tiferromagnetic layers the notions of AFM and FM lEC 
are not applicable, two types of co-linear correlated spin 
orientations in successive layers are still possible: iden- 
tical (in-phase) and reversed (out-of-phase), as shown in 
Fig. \m. Both types of lEC were observed in the experi- 
ment. 

All the (EuX),„/(PbX)„ SLs (where m and n denote 
the number of EuX and PbX layers, respectively) have 
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FIG. 2: Schematic view of two successive atomic layers of the 
EuS/PbS SL grown along [001] axis. The solid gray circles 
represent anions, the white circles cations - small circles are 
for ions in upper layer, the big ones for ions in the other. The 
distances rivjv from a cation to the four in-plane NN anions 
are shown by solid lines, by dotted lines the distances tnnn 
to the in-plane NNN cations 




FIG. 3: Schematic view of the crystallographic structure of 
the EuX/PbX (111) SLs. The ions in three successive cation 
planes are represented by gray, black and white circles, re- 
spectively. Only three NN anions, lying in the layer a\/3/6 
above the black cation plane, are shown by open squares. The 
distances to all twelve NNN are marked by dotted lines 

relatively small lattice mismatch. This mismatch, as well 
as the strains resulting from it, will be ignored in the fol- 
lowing, though the strains were shown to affect the mag- 
netic properties of the EuX layers, and, their transition 
temperatures to the paramagnetic state,LiE3. 

To discuss the spin coupling between the magnetic lay- 
ers one has to consider a SL magnetic elementary cell 
containing at least two such layers. In SLs grown along 
[001] crystallographic direction the situation is simple, 
as the stacking sequence is ABAB-type. Such stack- 
ing does not enlarge the size of the magnetic elemen- 
tary cell, whatever m + n value. In contrast, the stack- 
ing sequence ABCABC-type (compare Fig. ||) for both 
anions and cations, in SLs grown along [111] axis does 
enlarge the elementary cell when to -f n is not a multi- 
ple of 3. Thus, to limit the size of elementary cells, we 
consider only the (111) SLs with (n + m)/3 = integer. 



In the case of (001) SLs the primitive lattice vectors, 
which define our elementary cells are: ai = a/2[l, — 1, 0]; 
aa = a/2[l,l,0]; ag = a[0,0,m + 7i]. For (111) SLs, 
when the z-axis is chosen along the growth direction, we 
have ai = a/4[\/2, -^6, 0];a2 = a/4[V2, \/6, 0]; ag = 
(2a/-\/3)[0, 0, m + n] (the plane views of these cells are 
sketched in Fig. ^jand Fig. For both types of SL, our 
magnetic elementary cells contain, therefore, 2{m+n) an- 
ions, 2to magnetic cations and 2n nonmagnetic cations. 

In order to determine lEC in the above structures we 
compare the total valence-electron energies in two mag- 
netic SLs: with the in-phase and out-of-phase spin or- 
dering. For the n and m values typical for the exper- 
imentally studied SLs the elementary cells contain sev- 
eral tens of atoms. In view of this complexity, we decide 
to use the simplest calculation scheme still leading to 
fairly realistic band structure, namely, an empirical tight- 
binding method. Even though the one-electron methods 
are not designed for the total-energy calculations, the 
small spin-dependent changes in the total energy should 
be described adequately within this approach. 

III. TIGHT-BINDING MODEL 

To construct the empirical tight-binding Hamiltonian 
matrix one has to select the set of atomic orbitals for ev- 
ery type of involved ions and to specify the range of the 
ion-ion interactions. This selection is always a compro- 
mise between the best description of the band structure 
and the minimization of the Hamiltonian matrix dimen- 
sions and of the number of parameters used. In the fol- 
lowing we assume that the proper description of SL band 
structure is reached, when the Hamiltonian reproduces in 
the n = and m = limits the known band structures 
of the bulk constituent magnetic and nonmagnetic mate- 
rials, respectively. This criterion determines in principle 
the selection of the ionic orbitals and gives the values of 
the parameters, all but those characterizing the interac- 
tion between magnetic and nonmagnetic cations. 

To calculate the band structure of the lead chalco- 
genides (PbS and PbTe) we took into account s and p 
orbitals for both anions and cations, which lead to 8x8 
Hamiltonian matrix. We allowed for the s-s, s-p and p-p 
anion-cation nearest-neighbor (NN) interactions and the 
anion-anion and cation-cation p-p next-nearest neighbor 
(NNN) interactions. It turned out that the band struc- 
ture can be reproduced much better when we include 
also, by second-order perturbation, the interactions of 
p-orbitals with the three NN d-orbitals belonging to the 
F2 representation. The values of the parameters describ- 
ing all these interactions and the values of the on-site 
orbital energies were determined by a minimization 
procedure, in which the band structure was fitted to the 
energies inJJrie high symmetry points of the BZ, taken 
fronxj andll3. The obtained energy bands for PbS and 
PbTe along the symmetry axes of the BZ, are presented 
in Fig. ^. 
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FIG. 4: Model band structures of PbS and PbTe 



In the other limit, for europium chalcogenides (EuS 
and EuTe), to describe the cations we take explicitly 
the one s- and five d-orbitals, whereas the anions are 
described as before by s- and p-orbitals. The NN in- 
teraction involving the anion p-orbitals and cation s- 
and d-orbitals as well as the NNN cation-cation d-d and 
anion-anion p-p interaction were included in the 10x10 
Hamiltonian. The s-anion-s-cation, s-anion-d-cation in- 
teractions, turned out to be less important and were ne- 
glected. Instead, we included, again by second order per- 
turbation, the hybridization of anion p(|)-orbitals with 
the cation f(|)-orbitals - this was necessary for repro- 
ducing the spin splittings of the valence bands in the 
ferromagnetic EuS (we neglected the hybridization with 
the energetically distant f(J,) band). To reproduce the 
spin splittings in the EuS conduction bands the on-site 
exchange constants Js and Jd had to be introduced. The 
band structure of EuS, presented in Fig. ^ was obtained 
with the parameters fitted to the results of the APW 
spin-polarized calculations reported inta. The presented 
in Fig. H band structure of EuTe was obtained with the 
parameter values extrapolated from the values for EuS 
and EuSe by exploiting the chemical trends in europium 
chalcogenides. The elementary cell of the antiferromag- 
netic EuTe has a twice larger volume and completely dif- 
ferent shape than the one of the ferromagnetic EuS - to 
facilitate the comparison with EuS, we present the band 
structure of EuTe in the paramagnetic phase. 

In the above, the number of independent fitting pa- 
rameters was ipartially reduced according to the Har- 
rison relations^, e.g., instead of two NNN interatomic 



matrix elements ppa and pp'K, we used ppa = —Appn. 
It has to be also noted, that in all calculations we ne- 
glected the spin-orbit terms, known to be important in 
lead chalcogenides. These terms would increase the num- 
ber of model parameters and double the matrix dimen- 
sions and would, therefore, pose a problem in the case of 
the SLs. Fortunately, we are mainly interested in valence 
bands, for which the spin-orbit is much less important 
than for the conduction bands. In the magnetic elemen- 
tary cell of the (EuX)m/(PbX)ji SL there are 4(m -I- n) 
nonequivalent ions. Seven orbitals (s,p,d) have to be 
taken into account for each anion and each nonmagnetic 
cation and thirteen (s,d, and f) for each magnetic cation. 
In principle, the SL tight-binding Hamiltonian is there- 
fore a (40to -I- 28n)x(40TO + 28n) matrix, which has to be 
completed by the nonmagnetic cation-magnetic cation 
interactions. The constants describing these interaction 
can not be inferred from the m = and n = limits, 
they were estimated from Harrison's formula for inter- 
atomic matrix elements. The matrix was then perturba- 
tionally reduced to the (20m-|- 16n)x(20m-|- 16n) matrix. 
To determine the small difference between the large to- 
tal energies of the valence electrons in the two, in-phase 
and out-of-phase, spin configurations, we did not calcu- 
late these energies separately. Instead, the difference be- 
tween the two energies of the valence electrons was cal- 
culated at a given k-point, after the numerical diagonal- 
ization of the two corresponding Hamiltonian matrices 
and summation of the occupied states' energies. The 
results, obtained for a mesh of k-points (16000 points in 
the case of the tetragonal BZ of the (100) SLs, and 17280 
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in the case of the hexagonal BZ of the (111) SLs), were 
then used in the triple Simpson procedure for integrat- 
ing over the entire BZ. The calculations were performed 



for (EuS)„/(PbS)„ SLs grown in both [001] and [111] di- 
rections and for (EuTe),„/(PbTe)„ SLs grown along the 
[111] axis. 




We denote by AE the absolute value of the energy 
difference between the in-phase and out-of-phase spin 
configuration per unit surface of the layer. AE can 
be regarded as a measure of the strength of the inter- 
layer spin coupling in the SLs - for ferromagnetic struc.- 
tures it can be expressed in terms of the constant JiO, 
commonly used to characterize the lEC in metallic mag- 
netic/nonmagnetic trilayers by the relation: AE — 4| Ji| 
(here the factor 4, instead of 2, accounts for the fact, that 
in SL each magnetic layer is coupled to two neighboring 
layers) . 

The sign of the calculated energy difference determines 
the spin configuration in consecutive magnetic layers. In 
the ferromagnetic, EuS-based, SLs the out-of-phase spin 
configuration in consecutive magnetic layers is energet- 
ically preferred, so that lEC in these structures has an 
ANTIFERROMAGNETIC character, in agreement with 
the experiment. For the antiferromagnetic, EuTe-based 
SLs the situation is more complicated: for odd number 
m of spin planes in the magnetic layer the out-of-phase 
configuration has the lower energy, whereas for even m 
it is the in-phase configuration, which is energetically fa- 
vored. Thus, one can notice that in all studied SLs the 
valence electron mediated lEC prefers the spin configu- 



ration with the opposite directions of spins at the two 
interfaces bordering the spacer. 

Many various m and n values were selected to study the 
range of the interlayer coupling and the lEC dependence 
on the thickness of the magnetic layer. It turned out, 
that in all SLs for fixed spacer thickness n the strength 
of lEC is almost independent on the magnetic layer thick- 
ness TO. This seems to prove that in the considered here 
SLs, composed of two semiconductors with very differ- 
ent energy gaps, the valence electron mediated lEC is 
essentially a surface effect. 

The dependence of the strength of the interlayer cou- 
pling on the spacer thickness n for all three studied 
types of SLs is presented in Fig. ^. We recall that for 
the [111] growth direction, we calculated lEC only for 
(n + to)/3 — integer, so that for these SLs in the figure 
the points for different n values do not correspond to the 
same to. As can be seen in Fig. ^, the strength of the 
coupling in all three cases decreases with the spacer thick- 
ness approximately exponentially. The strongest and the 
least rapidly decreasing lEC was obtained in the case 
of the FM EuS/PbS (DDI) SLs. The comparison of the 
results obtained with the same set of model parameters 
for the two different types of EuS/PbS SLs, (i.e., the 
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FIG. 6: The interlayer exchange constant Ji (left Y-axis) as 
a function of the spacer thickness for ferromagnetic EuS/PbS 
(001) (squares) and (111) (triangles) SLs. For the antiferro- 
magnetic EuTe/PbTe SLs (circles) the absolute value of the 
energy difference AE was divided by 4 (at the right Y-axis) 
for comparison with the FM case 



grown along [001] and [111] axis) indicates that the va- 
lence electron mediated lEC depends strongly on the lat- 
tice geometry. For example, for the PbS spacer n = 2 the 
obtained coupling between EuS (111) magnetic layers is 
about five times weaker than that between EuS (001) 
layers. Moreover, in the (111) SLs the strength of the 
coupling decreases more quickly with n. The small regu- 
lar deviations from the smooth dependence of Ji vs. n, 
which can be best seen for every second n in the (001) 
case, but also for every third n in the (111) results, reflect 
the periodic effects of stacking. 

The lEC calculated for the AFM EuTe/PbTe (111) SL 
turned out to be stronger than that for FM EuS/PbS 
(111) SL (see Fig. ||). The difference in the band pa- 
rameters by itself does not explain this result - the cal- 
culations performed for AFM and FM (111) SLs with 
identical sets of band parameters have shown that the 
coupling between the AFM layers is approximately two 
times stronger. This indicates that the valence electron 
mediated lEC is not a purely surface effect. In view of the 
present results for lEC in antiferromagnetic EiiTe/PbTe 
SLs, we note that the simple ID model inliS overesti- 
mated the strength of lEC for large spacer thicknesses 
- now the decrease of lEC with n is much more rapid; 
the character of the interlayer coupling was however de- 
scribed properly. 



IV. COMPARISON WITH THE 
EXPERIMENTAL DATA 

A. Ferromagnetic EuS/PbS superlattices 

The lEC was observed in ferromagnetic EuS/PbS 
(001) SLSpjby magneticQ and neutron diffraction and 
reflectivityQ methods. For the samples with the thin 



i.e., 4.5A (1.5 monolayers, probably a 
2 and n — 1) and lOA (ca n = 3), 
the antiferromagnetic interlayer coupling was observed in 
the magnetic moment measurements; the magnetic meth- 
ods did not reveal any lEC in the samples with larger 
spacer thicknesses. In the neutronographic experiments 
the AFM lEC was confirmed in the above two samples, 
but it was also observed in the sample with 23A PbS 
layers. Further measurements in the external magnetic 
fields, parallel to the layers, allowed to estimate the ex- 
perimental strength of the coupling constant Ji , from the 
magaetic field B erasing the AFM neutron refiectivity 
peakEj. This was possible for the two samples with the 
thinner spacers, but not for the sample with the thicker 
spacer. For the latter, the field-induced changes in the 
neutronographic spectra were irreversible, what suggests 
that in this case the lEC was weaker than the magnetic 
anisotropics. The estimated experimental values of Ji 
are: 0.063, 0.031 and 0.019 (in mJ/m^), for n=l, 2 and 
3, respectively,E3. The corresponding theoretical values 
obtained from our model are: 0.77, 0.33, 0.18 mJ/m^. 
Thus, one can conclude, that the model of valence elec- 
tron mediated lEC describes properly the sign and the 
rate of the decrease of the coupling with the spacer thick- 
ness, but overestimates the strength of the coupling. The 
fact that the theoretical results obtained for crystallo- 
graphicaly perfect SLs lead to exchange constants order 
of magnitude larger than those observed for the real mul- 
tilayer structures, is probably due to the interface diffu- 
sion, which in the case of metallic structures wjas shown 
to reduce significantly the strength of the lECn. 



B. Antiferromagnetic EuTe/PbTe superlattices 

Unfortunately, for the AFM type of SLs there are no 
experiments, which provide direct information about the 
strength of the coupling. The evidences of the existence 
of the coupling between the AFM EuTe layers come from 
the satellite structure of the neutron diffraction spectra, 
seen in a variety of ElulCe/PbTe SLs consisting of several 
hundreds of periods ,ErE3. The detailed analysis of the 
shapes of the satellite lines in the neutronographic spec- 
tra indicates that in these SLs the EuTe layers are not 
entirely coupled, but only partially correlatecLthe less 
the thicker are the PbTe spacer layers of the SLl3. Under 
the strong assumption that the structures are morpholog- 
ically perfect, with the same m and n values throughout 
the entire SL, this degree of correlation can be quanti- 
tatively determined. Under the same assumption, the 
analysis of the satellites positions allows to distinguish 
which spin configuration, the in-phase or out-of-phase, is 
dominating. 

The observed spectra for the SLs with nominally even 
m and even n reveal the preference for the in-phase spin 
configurations, whereas those for the SLs with odd m 
and even n exhibit the preference for the out-of-phase 
configuration, both in agreement with the predictions of 
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the present model. For the case of even m and odd n 
there are no available data. Finally, for the samples with 
TO and n both odd, the neutron diffraction spectra seem 
to indicate that the in-phase configuration is preferred, 
contrary to the theoretical predictions. However, the in- 
phase spin configuration for SLs with odd number of spin 
planes in each antiferromagnetic layer should exhibit fcr- 
rimagnetic properties, i.e., lead to a significant net mag- 
netic moment of ferrimagnetic domains. No such mag- 
netic moments were detected in these samplesEZl. These 
somewhat confusing results seem to indicate that both 
chemical and magnetic structures of the studied SLs are 
not perfect enough. New technological and experimen- 
tal efforts to observe the lEC in EuTe/PbTe SLs with 
smaller number of SL periods, i.e., in SLs with better 
controlled m and n values, are undertaken. 

In conclusion, we have shown, within a 3D tight bind- 
ing model, that the valence electron mediated interlayer 
exchange coupling explains the AFM coupling between 
the FM layers observed in EuS/PbS (001) SLs with nar- 
row PbS spacers. The strength of the calculated coupling 



depends strongly on the lattice geometry and decreases 
approximately exponentially with the spacer thickness 
n. For a given type of SL it is almost independent on 
the number to of the spin planes within each magnetic 
layer. These features distinguish the considered mecha- 
nism from another mechanism of AFM coupling between 
the FM layers, namely the dipolar coupling possible in 
multilayer structures with tiny magnetic domainsEEI. 

The valence electron mediated interlayer coupling is, 
up to now, the only effective mechanism capable to ex- 
plain the origin of the interlayer correlations observed 
in the antiferromagnetic EuTe/PbTe SLs with no local- 
ized impurity states. The current, not complete under- 
standing of the experimental data for AFM SLs, does 
not allow, however, to draw definite conclusions about 
the comparison between the details of the experimental 
and theoretical results. 
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